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x-dimension of rectangular waveguide 
dominant mode reflection coefficient 

measured reflection coefficient 
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complex waveguide modal amplitude of mode 
y-dimension of rectangular waveguide 
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free-space wave number 
global finite element matrix 

element matrix for single tetrahedral element 

surface area over plane P l 

vector testing function 
dominant mode 

excitation column vector 

x-,y-, and z-dimensions of rectangular sample 

Cartesian Coordinate system 

modal admittance of p th mode 

permittivity of free-space 

permeability of free-space 

relative permittivity of medium in Region II 

relative permeability of medium in Region II 

real part of e r 

imaginary part of e r 
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small increment in £ " 
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real part of \x r 
imaginary part of |i . 
small quantities 

propagation constant for p th mode 
angular frequency 
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Abstract 


A simple waveguide measurement technique is presented to determine the 
complex dielectric constant of a dielectric material. The dielectric sample is loaded in a 
shorted x-band rectangular waveguide. Using a network analyzer, the reflection coefficient of 
the shorted waveguide (loaded with sample ) is measured. Using the Finite Element 
Method (FEM) the exact reflection coefficient of the shorted waveguide (loaded with the 
sample) is determined as a function of the dielectric constant. Matching the measured value 
of the reflection coefficient with the reflection coefficient calculated using the FEM utilizing 
the Newton-Raphson Method, an estimate of the dielectric constant of a dielectric material is 
obtained. A comparison of estimated values of dielectric constant obtained from simple 
waveguide modal theory and the present approach is presented. 

I Introduction 

Application of materials in aerospace, microwave, microelectronics and 
communication industries requires the exact knowledge of material parameters such as 
permittivity and permeability. Over the years many methods have been developed and used for 
measuring permittivity (e ', E ") and permeability ( (i ") of materials [1], The most accurate 

measurement at high frequency can be done using the high Q resonant cavity technique [2]. 
However, the main disadvantage of the cavity method is that the measured results are applicable 
only over a narrow frequency band [2]. Electric properties of material over wide range of 
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frequencies can be done with less accuracy using the transmission line methods. In the 
transmission line method an isotropic material sample with specific length is positioned in a 
transmission line and the (e e ") and (p ; .', pT') are determined from the measured reflection 

and transmission coefficients. The material sample used in these measurements is usually of a 
cross section which is the same as that of the trasmission line. The uniform cross section of the 
sample is selected so that a dominant mode analysis is sufficient and accurate for measuring the 
material constants. However, when the sample selected is not of uniform cross section or the 
sample occupies only a part of trasmission line cross section, then the complete modal analysis is 
required to measure accurately the material properties. The complete modal analysis is quite 
complicated, if not impossible. In such cases, when the sample cross section is different from that 
of the transmission line, a numerical method such as the FEM instead of the modal analysis is 
much easier to implement to obtain material properties. 

In the present report, the FEM is proposed to estimate complex permittivity of 
material using a terminated rectangular waveguide. The method described here may be easily 
extended to estimate the complex permeability of material. The material sample of specific 
length but of arbitrary cross section is assumed to be present in a shorted rectangular waveguide. 
The reflection coefficient at some arbitrary selected reference plane in the rectangular waveguide 
is measured at a given freqeuency. Since only determination of permittivity is required, a single 
reflection coefficent measurement suffices. The reflection coefficient at the given frequency is 
also calculated as a function of (e ', e ") using the FEM. From the calculated and measured 

values of reflection coefficients and the use of Newton-Raphson Method [3], the complex 
permittivity of the given material is then determined. The complex permittivities of Teflon and 
Plexiglass obtained using the present techniques are compared with the values obtained using the 
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standard software available with the ftp ( PCervCett Packard )- 8510 Network Analyze [6-7]. 


II Theory 


(a) Direct Problem 


In this section the FEM is used to determine the reflection coefficient of a short 
circuited rectangular transmission line loaded with an arbitrary shaped dielectric material. Figure 
1 shows a terminated rectangular waveguide with a dielectric sample of arbitrary cross section. It 
is assumed that the waveguide is excited by a dominant TE 10 mode from the right and the 
reflection coefficient is measured at the reference plane P| as shown in figure 1(a). For the 
purpose of analysis the problem is divided into two regions: Region I ( z < 0 ) and Region II 
(z > 0 ). Using the waveguide vector modal functions, the transverse electromagnetic fields in the 
region I are expressed as [4] 
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In deriving equation (1) and (2) it is assumed that only the dominant mode is incident on the inter- 

face P] and a p is the amplitude of reflected modes at z=0 plane. e 0 (x, y) and ho (x, y) are the 

— ^ ^ 

rectangular waveguide vector modal functions for the dominant mode and e p (x,y) h (x,y) are 
the rectangular waveguide vector modal function for p th mode [4]. Y p and appearing in 
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th 

equations (1) and (2) are the characteristic admittance and propagation constant for p mode and 
are given by 



where m, n are the mode indices, a, b are the x- and y-dimensions, respectively, of the 
rectangular waveguide, and k Q is the free-space wave number. The unknown complex modal 

amplitude a p may be obtained in terms of the transverse electric field over the plane P 1 as 
follows 
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n 


-> 

E 


over S, 


-> 

E 


over Sj 


• e 0 ds 


• e p ds 


(3) 

(4) 


— ^ 

where E 


over s is a tangential electric field over the surface S l at the reference plane P { . 

The electromagnetic field inside Region II is obtained using the FEM formulation [5]. 


The vector wave equation for the E field is given by 


, 1 ^ n \ ( 2 V" 

Vx - ■ VxE - k 0 e )E = 0 




(5) 


/V /V 

Multiplying equation (5) by a testing function T (x, y, z) = xT x + yT y + 'z,T_ and integrating over 
the volume of the Region II, equation (5) yields with simple vector manipulations 
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where V is the volume of Region II and h is the unit normal vector to the surface S l drawn out- 
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ward with respect to the volume V . Using VxE = -ja\x 0 H , equation (6) may be rewritten as 
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In order to solve the equation (9), the volume enclosed by Region II is discretized by 


using first-order tetrahedral elements such as the one shown in figure 2. The electric field in a sin- 



gle tetrahedral element is represented as 

-»// -> 
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where Z? , m=l,2,3..6 are the six complex amplitudes of electric field associated with the six 

— ^ 

edges of the tetrahedron, and W m (x, y, z) is the vector basis function associated with the mth 

— ^ 

edge of the tetrahedron. Detail derivation for the expression for W m (x, y, z) is given in refer- 
ence [5]. Substituting equation (10) into equation (9), integration over the volume of one tetrahe- 
dron results in 
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The above equation can be written in a matrix form 
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These element matrices can be assembled over all the tetrahedron elements in Region II to obtain 


a global matrix equation 



The solution vector of the matrix equation (15) is then used to determine reflection 
coefficient at the reference plane P l as 
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(b) Rectangular Waveguide Measurement System 


(16) 


The reflection coefficient a 0 ' due to a terminated reactangular waveguide loaded with 

a given material sample piece can be measured using the procedure described in [6]. Assuming 
the sample piece occupies the entire cross section of the waveguide, an algorithm which uses the 
Nicholson-Ross Technique [6-7] is used to determine the complex permittivity of the sample. 
However, the algorithm which is based on the Nicholson-Ross Technique cannot be used when 
the sample piece occupies part of the cross section of rectangular waveguide. Also the 
Nicholson-Ross Technique can only be used for a through measurement and not for a 
measurement with terminated waveguide. When the dielectric sample is of arbitrary shape the 
procedure described in the next section is used. 

(c) Inverse Problem 

This section presents computations of complex dielectric constant of a given sample 
piece from a one port measurement of the reflection coefficient. From the given geometry of the 
sample and its position in the terminated rectangular waveguide, the reflection coefficient 
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a Q (e r ', e") is calculated using the FEM for assumed values of (£ r \ e") . If a Q ' is the measured 
reflection coefficient then the error in the calculated value of reflection coefficient is 
a Q (&;, £") - a Q ' . Writing the error in real and imaginary part we get 

A (e;, = real (a Q (£/, £/') - a 0 ') (17) 

f 2 (e;, e r ") = imag (a Q (£/, £ .") - a 0 ') (18) 

If (£/,£/') are incremented by small values to (e; + de;, e” r + de" r ) such that 
A (£ r ' + de;, e" r + de" r ) and f 2 (e; + <r/£ ; .', £" r + de" r ) are simultaneously zero then we can write 
following matrix equation [3] 
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From the solution of equation (18), new modified values of (£ r ', £ ; ") are obtained as 

(£/) = £/ + de; ( 20 ) 

(£") =£" + j£" (21) 

v r > new r r v J 

The reflection coefficient a 0 (£ r ', £ r ") with modified values of (£ r \ £ r ") is again calculated 
using the FEM procedure. With the new value of a Q (£ ; \ e") computation through equations 
(17)- (19) are repeated. The above procedure is repeated until required convergence is obtained 
(i.e. <r/£ ; .'<A 1 and de" < A 2 , where Aj and A 7 are preselected small quantities ). The 

procedure described above will converge to the true value of complex permittivity if the first 


11 



choice of (£ ; .', £ ") is close to the true values of complex permittivity. 

Ill Numerical Results 

For numerical result , samples made from Teflon and Plexiglass materials are 
considered. First, samples which occupy the entire volume enclosed by Region II are considered. 
This is done for the purpose of comparison with results obtained by the present technique and the 
results obtained by standard software available with the ftp (- 8510 Network Analyzer. Numerical 
results obtained from samples which do not occupy entire volume of the Region II are also 
presented. 

First, a sample piece of rectangular shape with dimension 
w x = 2.29 cm, w = 1.0 cm, w_ = 0.947 cm was cut from a Teflon sheet. This size was 
selected so that the sample piece occupies entire volume of the Region II. The reflection 
coefficient at reference plane P l was measured by placing the sample piece in a x-band 

rectangular waveguide and using ftp - 8510 network analyzer over the frequency band 8.2GHz- 
12.40GHz. Using the procedure described in section II, the complex permittivity of the Teflon 
sample is calculated and shown in figure 3. For the same sample complex permittivity using the 
Nicholson-Ross Technique is also calculated and presented in figure 3. The two results agree 
well except at a frequency close to 12 GHz where small resonance is detected by the FEM 
procedure. Because of very low loss characteristic of the Teflon material, the estimation of 
imaginary part of complex permittivity is not very reliable. The complex permittivity of 
Plexiglass is also obatined using the present procedure and presented in figure 4 along with the 
results obtained by Nicholson-Ross Technique [7]. The two results are in good agreement with 
each other. 
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In order to validate the present method for a material sample for an arbitrary shaped 


sample, a Teflon sample with dimension w x = 1.58 cm, w = 0.785cm, w„ = 0.632cm is cut 

from a Teflon sheet. The sample is then placed in a rectangular waveguide for measurement of 
reflection coefficient as shown in figure 5. From the measured value of reflection coefficient and 
following the procedure described in section II the complex permittivity is calculated and 
presented in figure 6. For comparison, the complex permittivity obtained using full size sample 
(i.e. w x = 2.29 cm, w = 1.0 cm, w_ = 0.947 cm ) is also presented in figure 6. From figure 6 

it may be concluded that the FEM procedure can be used to determine the complex permittivity of 
dielectric material using an arbitrarily shaped sample. 

Figure 7 shows the real part of complex permittivity of Plexiglass obtained using a 
sample with dimension w x = 1.58 cm, w y = 0.785cm, w. = 0.632cm and placed in Region II 

as shown in figure 5. The real part of complex permittivity of Plexiglass obtained using full size 
sample is also presented in figure 7. The two results are in good agreement within a range of 
uncertainty specified by manufacturer. 

IV Conclusion 

A FEM procedure in conjunction with the Newton-Raphson Method has been 
presented to determine complex permittivity of a dielectric material using an arbitrarily shaped 
sample. The arbitrarily shaped sample of a given dielectric material is placed in a terminated x- 
band rectangular waveguide. The reflection coefficent at a reference plane is measured with h ip- 
8510 Network Analyzer. For the same configuration of a terminated x-band waveguide loaded 
with the sample piece, the reflection coefficient is calculated as a function of complex dielectric 
constant using the FEM . The Newton-Raphson Method is then used to determine the complex 
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dielectric constant by matching the calculated value with the measured value. The measured 

values of complex permittivity of Teflon and Plexiglass using the FEM method are in good 

agreement with the results obtained by the Nicholson-Ross Technique. 
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Figure 1(a) Longitudinal view of rectangular waveguide 
loaded with dielectric sample 



Figure 1(b) Cross-sectional view of rectangular 
waveguide 


Figure 1 Geometry of rectangular waveguide excited by TE 10 
mode 
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Figure 2 First order tetrahedra element 
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Figure 3(b) Imaginary part of complex permittivity of a Teflon sample over 
the x-band calculated using two different methods. 



Figure 4(a) Real part of complex permittivity of a plexiglass sample over the 
x-band calculated using two different methods. 
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Figure 4(b) Imaginary part of complex permittivity of a plexiglass sample over 
the x-band calculated using two different methods. 
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Figure 5 Geometry of the x-band rectangular waveguide holding dielectric 
sample piece. 
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Figure 6 Real part of complex permittivity of a Teflon sample obtained 
using undersize sample. 
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